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This work studies jet-like electromagnetic configurations surrounding a slowly-spinning black-hole 
immersed in a uniformly magnetized force-free plasma. In the first part of this work, we present 
a family of stationary solutions that are jet-capable. While these solutions all satisfy the force- 
free equations and the appropriate boundary conditions, our numerical experiments show a unique 
relaxed state starting from different initial data, and so one member of the family is likely preferred 
over the others. In the second part of this work, we analyze the perturbations of this family of 
jet-like solutions, and show that the perturbative modes exhibit a similar split into the trapped and 
traveling categories previously found for perturbed Blandford-Znajek solutions. In the eikonal limit, 
the trapped modes can be identihed with the fast magnetosonic waves in the force-free plasma and 
the traveling waves are essentially the Alfven waves. Moreover, within the scope of our analysis, we 
have not seen signs of unstable modes at the light-crossing timescale of the system, within which 
the numerical relaxation process occurs. This observation disfavors mode instability as the selection 
mechanism for picking out a preferred solution. Consequently, our analytical study is unable to 
definitively select a particular solution out of the family to serve as the aforementioned preferred 
hnal state. This remains an interesting open problem. 

PACS numbers: 04.70.Bw, 94.30.cq, 46.15.Ff 


I. INTRODUCTION 


The potential role of magnetospheres around compact 
objects in helping power energetic phenomena has long 
been recognized m- Examples of such phenomena in¬ 
clude pulsars and active galactic nuclei (AGN) which rep¬ 
resent exciting laboratories to understand physics in ex¬ 
treme regimes. To date, significant insights into phenom¬ 
ena tied to magnetospheres have been gained through 
combined theoretical (analytical and numerical) efforts 
and contrasting with observations. On the theoretical 
fronts, considerable efforts have concentrated on the de¬ 
tails of the magnetosphere dynamic as it interacts with 
compact objects by making use of the force-free approx¬ 
imation (see, e.g. [4ll24jl. Of particular relevance to 
AGNs -as well as other related phenomena like gamma 
ray bursts, ultra-luminous X-ray binaries, etc-, is the 
question of how black holes power jets. Key insights 
towards answering this question have been provided by 
complex simulations as the inherent complexity of the 
relevant equations of motion, together with particularly 
extreme physical set-up (which requires dealing with a 
black hole), hinders analytical work. Beyond particular 
details, the common message from simulations is that 
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the magnetosphere can tap into the black hole’s kinetic 
energy reservoir and induce a strong Poynting flux from 
the magnetosphere. This understanding, complemented 
with further studies as to the stability of jets, interaction 
with a possible accretion disk, etc. has been instrumental 
in putting forward particular models of how jets may be 
launched (see e.g. [5SJ[5S] and references cited therein). 

Despite the aforementioned advances, it is arguably 
desirable to also gain access to black hole-magnetosphere 
interaction through analytical means. This would allow a 
more convenient way to explore questions like: What col¬ 
limates the emitted radiation? What is the electromag¬ 
netic field conhguration within the jet? How do the elec¬ 
tromagnetic and plasma perturbations propagate within 
the jets? These questions are ultimately intimatedly tied 
to observations. 

To help answer these questions, we here carry out ana¬ 
lytical studies complemented by numerical investigations 
within the force-free approximation and study the cor¬ 
responding magnetosphere as it interacts with a (slowly) 
spinning black hole. That such approximation is relevant 
in this context was pointed out over four decades ago 
in[T]. At the core of this observation is the fact that a 
particle-production cascade will take place and the mag¬ 
netosphere will be composed of a tenuous plasma which 
will reach equilibrium when the electric field becomes or¬ 
thogonal to the magnetic field. The electromagnetic en¬ 
ergy density of the plasma dominates the plasma mass 
density [3], and it is expected to remain so over time; 
correspondingly the transfer of energy and momentum 
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into the plasma is negligible. In mathematical terms, 

Ff,.f=0, ( 1 . 1 ) 

where F^i, is the electromagnetic field tensor and j^, is 
the 4-current density. This is the so called force-free con¬ 
dition maul] which results from = 0. 

In the present work, we focus particularly on revisiting 
the question of jets induced by spinning black holes. To 
simplify the analysis, we shall assume that the black hole 
region is threaded by an electromagnetic field configura¬ 
tion which, asymptotically, is purely magnetic and paral¬ 
lel to the black hole’s angular momentum. Consequently, 
the system has a rotational symmetry around the spin 
axis. In addition, we restrict ourselves to slowly-rotating 
black holes, and solve the force-free equations through 
a perturbative expansion (considerably generalizing the 
solutions found in laiZH!). Despite these approximations, 
we expect the physics discussed here to carry through to 
more general cases, for example with misaligned spins 
[29j or more rapidly spinning black holes, and expec- 
tectations already supported by numerical simulations 
(e.g. [IS1I2H]). We note that while prior analytical ef¬ 
forts have also concentrated in our case of interest [25] . 
we here adopt the powerful formalism of Euler potentials 
and the techniques of exterior calculus as discussed in 
Refs. |5nH55] and so can explore the possible solutions 
in broader terms. Our analysis will reveal the existence 
of a family of solutions describing jets which emanate 
from the black hole and propagate towards infinity along 
the symmetry axis. These jets are enclosed by a tube¬ 
like surface, within which energy and currents flow to 
infinity, and the electromagnetic field does not asymp¬ 
tote to the external constant field set up initially. Out¬ 
side of this tube, there is no charge, current or energy 
flux, and the field matches the external driver at infin¬ 
ity. An important observation is that all solutions within 
this family exhibit the same qualitative behavior we have 
just described, and we have not found a definitive selec¬ 
tion criteria as to which particular one is more physically 
realistic. In particular, none of these solutions appear 
to have an instability over the light-crossing timescale 
(across the tube), as suggested by the numerical investi¬ 
gation in Sec. m However, numerical studies of the jet 
solutions obtained with different initial data does sup¬ 
port the existence of a single relaxed state determined 
only by the black hole and the external field. Thus, it 
appears that such a criteria may very well exist and we 
are failing to impose it when solving for our stationary 
jet-like solutions. From our discussion, it will be clear 
that the task of selecting a particular solution is tightly 
connected with the resulting angular velocity of magnetic 
field lines and the development of a current sheet -where 
the force-free equations are inconsistent without remov¬ 
ing some amount of energy-. These are rather subtle and 
long-standing issues in the search for FFE solutions. We 
expect any successful future resolution of them to also 
provide important clarity to FFE dynamics in general. 


Having obtained such a family, we then investigate its 
perturbative mode structure. We identify a clean and 
simple separation into “trapped modes” and “traveling 
waves”, which are analogous to the findings made on the 
perturbed Blandford-Znajek configuration m- Interest¬ 
ingly, the trapped modes behave similarly to the vac¬ 
uum electromagnetic quasinormal modes (QNM) of black 
holes |55IH5| . In fact, in the short wavelength (eikonal) 
limit, these trapped modes become the fast-magnetosonic 
waves of the force-free plasma and satisfy exactly the 
same equations that the vacuum QNMs obey in the same 
limit. On the other hand, the traveling waves generically 
carry charge, and propagate inward towards the black 
hole or outward towards spatial infinity along the back¬ 
ground magnetic held lines. In the eikonal limit, they 
become the Alfven waves of the force-free plasma. These 
observations are consistent with the analysis of Uchida 
[82] . where he hnds fast magnetosonic waves and Alfven 
waves by examining local dispersion relations. Moreover, 
our study reveals no mode instability at the light-crossing 
timescale for any of these jet-like solutions, in contrast 
with the fast-variation that we observe numerically in 
Sec. |HI[ This suggests that stability can not be the cri¬ 
teria to single out the unique relaxed state observed in 
numerical simulations. 

We organize the paper as follows. In Sec. |H| we in¬ 
troduce the family of jet-like solutions and discuss their 
validity. In Sec. |HI[ we perform our numerical studies 
with different initial data. The hnal relaxed states are 
extracted and compared with each other. In Sec. |IV| we 
apply the formalism developed in Ref. m to analyze 
perturbations of these jet-like configurations, where dis¬ 
cussions on mode instability are made. Finally, we con¬ 
clude in Sec. jV] Throughout this manuscript we adopt 
the geometric units, setting the gravitational constant G 
and the speed of light c to one. 


II. STATIONARY JET SOLUTIONS 

In order to solve for the stationary held conhguration 
in the presence of a slowly-spinning black hole, we start 
with a background conhguration (the “Wald-type solu¬ 
tion”) for a Schwarzschild black hole, and then compute 
the changes required after a spin is introduced. We shall 
present a family of solutions which satisfy the force-free 
equations and the boundary conditions. 


A. The Wald solution for Schwarzschild black holes 

The Wald solution [44] describes the electromagnetic 
held distribution surrounding a axis-symmetric black 
hole immersed in a constant magnetic held. In our case, 
we start with the Wald conhguration for a Schwarzschild 
black hole. For simplicity, let us assume that the mag¬ 
netic held is oriented along the 2 direction (while a;, y 
directions are orthogonal to z), with its strength being 
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Bq. More specifically, The magnetic field in the Wald 
solution is (in orthonormal coordinate) 


B = Bo (cosde^ —-\/7sin0eg), (2.1) 

with / = 1 — 2M/r. The corresponding field tensor is 
F =B^dr Add — Bgdf A dcj) + BfdO A d(j) 

= ^d{rsm9)^ A d(j) 

=Bod{r sin 6 cos (j)) A d{r sin 6 sin (j)) 

=Bodx A dy . (2-2) 

In the language of Refs. [HOHSS] . the field tensor can 
be expressed in terms of the differentials of two “Euler 
potentials”, in the form of F = d(j)i A d(j) 2 , and these 
potentials are not unique. With the Wald solution, we 
can either choose the gauge in which the Euler potentials 
are ipi= r sin 9 cos (p and 1^2 = r sin 9 sin <j), or 



(/)i = sin^ 6*/2 , (f)2 = (j), (2.3) 

together with a multiplicative constant Bq. We will 
choose the latter in our calculations in the next sec¬ 
tion. Recall that the Wald solution was originally de¬ 
rived for vacuum (without plasma) Maxwell equations 
but in the case of a non-spinning black hole it satis¬ 
fies all the force-free requirements as there is no current, 
charge, or electric field in the spacetime. Therefore in 
the non-spinning case, it also provides a solution for a 
Schwarzschild black hole immersed in constantly magne¬ 
tized force-free plasma. 


FIG. 1: A schematic illustration for a jet-like solution 
launched from a rotating black hole. The magnet field lines 
originated at the 6 = 7r/2 ring of the horizon separate inner 
and outer regions of the jet. In the inner region, current and 
the angular velocity of the field lines are related by the Znajek 
boundary condition, whereas they are less constrained in the 
outer region. The red arrow labels the rotation direction of 
the black hole. 


given by Eq. (2.31 in the non-rotational limit, we shall 
have the following expansions 


B. Small spin expansion and jet solutions 

As shown in Ref. [3^, the force-free Euler potentials 
of a stationary, axis-symmetric spacetime can only take 
on a very restrictive form 

4>i ='4){r,9), (j )2 = i’2ir,9)+(j)-nF{'<P)t, (2.4) 


Ip = sin^ 9/2 + aip^^\r, 9) + 'ip^'^\r, 9) + 0{a ^), 

1/2 = a'tp!2^\r,9) + O(a ^), 

Qf =an^J:\ip)+0(a^). ( 2 . 6 ) 

We expand tp to higher orders in a here for reasons which 
will become clear later. In addition, the polar current is 
defined as 


where pj corresponds to the magnetic flux enclosed by the 
circle at constant r and 9, and fli? is the angular velocity 
of the held lines passing through that circle. In addition, 
the force free condition Eq. (I.I) turns into the pair of 
equations: 


d(pi A d* F = 0 , d(p2 Ad* F = 0 . 


(2.5) 


The potentials (2.3) are obviously in a form consistent 


with Eq. (2.4), and we will use them as the background 


conhguration. We can then turn on the black hole spin 
and investigate the change of Held distributions. As the 
spin of black hole a is assumed to be small, the deviation 
of the Euler potentials from their background values may 
be solved in a series expansion of a. Since (pi and (p 2 are 


* {dp> A dip2) = T^dt A d(p, 
2'k 


(2.7) 


where the operation with Hodge star * in a slowly- 
spinning Kerr spacetime is discussed in the Appendix 
of Ref. |37j- Using the Hrst expression of Eq. (2.5), 


this polar current can be shown to be a function of p) 
only [33], as the current Hows along the magnetic Held 
lines in the frame that co-rotate with the Held lines. 
Moreover, the polar current can also be expanded in a: 
I = al^^'^{pj) + O(a^), and we have 

a * [d(r^ sin^ 9) A dp}^'’] = ^—dt A d(p , (2.8) 

27r 
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and so 


d^|; A d'!/'2 = — * *{dijj A 

1 


27r / sin( 


d9 A dr. 


(2.9) 


The polar current I and the angular velocity ftp are 
not independent of each other - they are related by the 
Znajek boundary condition |45H47) . which requires that 
the held on the horizon be regular. In the current con¬ 
text, Znajek boundary condition reads 




= 2n{aQp'’ — sin0 


r+ 


r+ 


27r(afi^^ — VtH)p'^ 


r\ - 




( 2 . 10 ) 


where we have dehned p = r sin 0 as the cylindrical radial 
coordinate, and r+ = M + Vis the outer hori¬ 
zon radius and VIh = a/{2Mr+) is the horizon angular 
velocity. As and 12^^ are both function oiip Ki p^/2, 
they must be approximately functions of p as well. Con¬ 
sequently their functional relationship in the entire re¬ 
gion where p < r+ is the same as their relationship on 
the horizon. It is then clear that the Znajek boundary 
condition has naturally divided the spacetime into two 
parts (see Fig. for a schematic illustration). The in¬ 
ner part is enclosed by a tube, where the polar current 
and the field-line angular velocity is related with each 
other by Eq. (2.101. Outside of the tube these two quan¬ 
tities are less constrained. In principle, any solution of 
the form prescribed by Eq. |2.4| that satisfies the stream 
equation (i.e. the second force-free equation in 2.5 which 
we haven’t used. See next section) is admissible. 

There are arguments favoring zero ^Ip and polar cur¬ 
rent density in the outer region: In the inner region, 
the rotation of the held lines and the currents are both 
sourced by the rotation of the black hole, whereas in the 
outer region they have to be driven by sources at inhnity. 
As it is not physically motivated to have such sources at 
inhnity, we should expect flp = 0 and 1 = 0. In fact, as 
shown in 1361. we have 


where P identihes some curve in the (r, 9) plane. Physi¬ 
cally, the equations above imply that non-vanishing Qp 
and I would always introduce energy and angular mo¬ 
mentum huxes, which are injected (along the held lines) 
by the sources at inhnity and eventually escape out. As 
it is physically more nature to imposes outgoing bound¬ 
ary conditions at inhnity, such hux injections should not 
be allowed. While these arguments, strickly speaking, 
only apply to regions near the outer boundary, instead of 
the entire space outside of the tube, previous numerical 
simulations do favor zero Qp and constant I outside the 
jet-tube Hg. Therefore in the analysis below we mainly 
focus on the solution inside the tube. 


We note however, that there is one subtlety when we al¬ 
low the presence of a current sheet, where currents can be 
launched at, and move away from the equatorial plane. 
This may provide an additional driving mechanism for 
n p and / outside of the black-hole-driven tube, and effec¬ 
tively enlarge the jet-tube. This is likely what we observe 
in Fig. 1^ where the tube radius of the final relaxed state 
is clearly greater than r_|_. As the analytical description of 
the dynamical effect of a current sheet feeding back onto 
the force-free plasma remains elusive, we are not able to 
further constrain these current-sheet-driven regions. 

Finally, we note that although the discussion in this 
section relies on an expansion in a, we expect this basic 
picture of a tube dividing the magnetosphere into a jet 
region and an outside region to remain valid when a is not 
a small quantity. In particular, inside the jet, the Znajek 
boundary condition (without taking the small a limit as 
in Eq. 2.10) should still provide us with a relationship 
between / and flp. 


C. The stream equation 


The stream equation is derived from the force-free Con¬ 
or more precisely the second expression 
), and it is also often referred to as the 
, 115] . In a covariant lan- 


ditions (Eq. 1.1 
Eq. 


Grad-Shrafranov equation [151 SS] • 
guage it can be written as [55] 


II 


,-p 




= 0 , ( 2 . 12 ) 


where 


rj = d(j) — flp{ip)dt. (2-13) 

Given / and Up’s dependence on ip, the stream equa¬ 
tion determines the distribution of tp in the spacetime. 
In the case of a slowly-spinning Kerr background, it has 
been shown in Ref. [3] ^ that is zero and is 
sourced by and I^^\ In particular, for any Hp'^ 
and I^^\ as long as they satisfy the Znajek boundary 
condition, the resulting from the stream equation is 
well defined (bounded). Therefore the small-a-expansion 
analysis implies that there is a family of force-free elec¬ 
trodynamic solutions in the jet-tube, corresponding to 
different choices of the angular velocity function Up. As 
mentioned earlier, inside the jet-tube, the magnetic field 
does not asymptote to the external driving held, which 
is a fundamental difference between force-free electrody¬ 
namics and vacuum electromagnetism. In Sec. Ill we 
shall explore this observation numerically, but Hrst, let 
us discuss two physically interesting solutions. 


^ With a slightly different set of notations. We have the following 
identifications between the quantities of that paper and those 
utilized in the current presentation: X{r,9) <-> P^/2, x{r,9) o 
W{r,9) -H. n^p\ and Y(r,6») o lO)/l^2n). 
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D. Physically interesting solutions 


For this case, the field tensor inside the tube is 


We have now obtained a family of regular (on the hori¬ 
zon and elsewhere) solutions, which is too broad for the 
purpose of understanding jet physics (also see [SD] for 
discussions on different stationary solutions). Ideally, we 
would like to single out a unique physically realistic so¬ 
lution. However, as recent experiences with analytical 
searches for solutions in other spacetimes (such as in the 
near-extremal Kerr back holes [5T1 - I55] 1 show, being un¬ 
able to pick a unique solution is a natural state of affairs. 
Routinely, the condition of regularity is insufficient to pin 
down a unique solution. 

We mention in passing that a jet solution has been 
found in Refs. |54[ 155) . and their conclusion would seem 
to imply the solution found is unique. Motivated from 
the I vs. ilp relationship of the monopole solution, the 
jet solution obtained in these works is indeed valid and 
belongs to the family of solutions we have presented here. 
However, the additional constraint employed for singling 
out a unique solution in |54[l55j is motivated by the desire 
to link to the split monopole solution which need not be 
the physically natural condition. 

Despite the lack of any solid conditions to single out 
particular solutions from the family, we nevertheless note 
that there are solutions within this family possessing in¬ 
teresting physical properties. For instance, demanding 
that the solution maximizes the energy radiated d£/dt 
from the jet, it is clear from Eq. (2.11) and the Znajek 


boundary condition Eq. (2.10) that the optimal choice is 


flp = Dj:// 2. In this case, the field tensor inside the flux 
tube is given, to 0(a), by: 


F = 


ar'^ sin 


16/M3 
ar^ sin 0 cos 0 


dr A d0 


, ar sin^ 0 , 

dt A d0 -- 7 ^—dr A dt 


8M2 8M2 

sin 0 cos 0d0 A d(j) + r sin^ 0dr A dcj) (2.14) 


while outside the tube, we have the background “Wald 
solution” of Eq. (2.2). 


If one instead requires that the energy stored within 
the tube is minimized, we realize that only the part of 
the tube in the asymptotically fiat region needs to be 
considered, which accounts for most of the volume and 
contains most of the energy. In other words, we can 
simply require 


Pe cc = i?| 




1 -I- 


jBr 

(2^)2p2 






1 H-~ ~ p^) + 


(2.15) 


to be minimized, giving 


“ 2r2 -p2^^' 


(2.16) 


F = ^^v/4M2 -p2dr Ade 

sin 0 cos 0d0 Ad4> + r sin^ 0dr A dcj ), 

+r cos Oxdt A d0 — sin 0xdr A dt (2-17) 

where 

^ ^ 8M2 - p2 (2.18) 

III. NUMERICAL SEARCHES FOR THE 
RELAXED STATE 

A powerful way to analyze the steady state solutions 
possibly achieved by a system described by the force- 
free equations is to study numerical solutions obtained 
under different initial configurations and examine their 
late-time properties. To this end, we employ a numerical 
code described in [16] which has been thoroughly tested 
and employed to study different systems in e.g [5M55] . 
This code implements the force-free equations and en¬ 
sures sufficient resolution is achieved efficiently through 
the use of adaptive mesh refinement (AMR) via the HAD 
computational infrastructure that provides distributed, 
Berger-Oliger style AMR [Ml ISO] with full sub-cycling in 
time, together with an improved treatment of artificial 
boundaries |6T] . We also note that a Cartesian grid stuc- 
ture is employed without imposing any symmetries. As 
a result, truncation errors introduce the typical m = 4 
perturbation. As we will illustrate next, a rather unique 
steady state solution is achieved irrespective of the ini¬ 
tial configuration, provided the asymptotic conditions on 
the electromangetic field are the same. This observation, 
together with the fact that the m = 4 perturbation does 
not seem to affect the observed behavior reinforces the 
belief that it is attractive and stable. 

For simplicity, and since the electromagnetic field’s 
back-reaction on the spacetime is negligible, we adopt 
for the geometry a spinning black hole solution in terms 
of ingoing Eddington-Finkelstein coordinates. For the 
initial electromagnetic fields, we adopt the two physical 
cases described in section ITl Dl which we label case PI and 
case P2 for the radiated energy maximizing and minimiz¬ 
ing energy stored configurations respectively. As well, we 
use an additional initial data with Ea = Fabn^ = 0 and 
Ba = *Fabn^ = Bq51 (case Al). 

We note that PI, P2 of these initial configurations are 
stationary solutions of the force-free equations, while Al 
is not. Nevertheless, in all cases, one observes an Alfven- 
wave-like transient that is radiated away from the equa¬ 
torial plane and then, after a few crossing times, the so¬ 
lution is seen in all cases to relax to the same solution. 
As an illustration, Fig.j^shows the value of Vlp obtained 
after the relaxed state is reached for the PI and P2 cases, 
together with a higher resolution run for case Al (Fig. 
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EH 


FIG. 2: Qf versus x at z = 7M at different early times in a 
spinning black hole with a/M — 0.1 for cases PI, and P2. 



FIG. 3: versus x at z = 7M at different times in a spinning 

black hole with a/M = 0.1 for cases PI, and P2 and also the 
case A1 at / = 120M. Despite the initial differences, all three 
solutions relax to a single one. 



FIG. 4: L 2 norm of the difference between the vaues of Df 
obtained with the two “physical” cases for a spinning black 
hole with a/M — 0.1. 


split-monopole configuration, and here we apply similar 
techniques to the axis-symmetric and stationary jet-like 
solutions discussed in Sec. [l^ We shall first explain the 
perturbation method, and demonstrate its usage by find¬ 
ing the modes of general axis-symmetric and stationary 
solutions in flat spacetime. After that we shall tackle the 
more difficult task of analyzing the mode structure of jet¬ 
like solutions of the slowly-rotating Kerr spacetime. At 
the end of the section, we will comment on the stability 
of this family of solutions. 


A. Formulation 


For any stationary and axis-symmetric solution de¬ 
scribed by Eq. (2.4), the perturbed force-free configu¬ 
rations can be written as 


shows the earlier evolution of this quantity). Clearly, 
despite different initial configurations adopted, there is 
an agreement in the relaxed states achieved by all cases. 
Moreover, the final solution is axisymmetric as expected, 
which is illustrated in Fig. |^that plots the (norm of the) 
difference between the 12 p values on the a; = 0 and y = 0 
planes, as a function of time and for both the PI and P2 
cases. As time progresses, the differences decay exponen¬ 
tially to zero. 


IV. THE MODE STRUCTURE OF THE JET 
SOLUTIONS 

In addition to the stationary configurations, it is also 
interesting to study the dynamical evolution of black hole 
magnetospheres and the jets. In fact, many radio obser¬ 
vations of quasars and pulsars have revealed significant 
time-dependent variations (see e.g. Refs. [621 - I65] '). which 
invite further studies of the generation and propagation 
mechanisms of perturbations within the jets. In Ref. m 
we studied the mode structure of the Blandford-Znajek 


F = Fp + 6F = + eda) A {d(j )2 + ed/3), 

where a, /3 are the perturbative Euler potentials and e is 
a flag that helps us track the order of these perturbative 
fields. The perturbative part of the field tensor is given 

by 


SF = da A d(j )2 + d</'i A d/3 . 


At 0{e), the force-free Eqs. (2.5) then become 


da Ad* F + d(j>i Ad* SF = 0 , 
d/3 A d * F + d(/)2 Ad* SF = 0 . 


(4.1) 


(4.2) 

(4.3) 


Before we attempt to solve these equations, we note 
that when we substitute in the background Euler poten¬ 
tial (/)2 (2.5) for the slowly rotating solutions, the force- 
free equations above pick up an explicit t dependence 
tied to our gauge choice. This will of course not feed 
into the final results for the perturbed field tensor, which 
should respect the time-translational symmetry of the 
background solution, but is certainly undesirable. It is 
both preferable and convenient to refine the perturbing 
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field /?, so that the resulting force-free field tensor be¬ 
comes 


F = Fb + 5F 

= {dij) da) A {d\j }2 — d{VLFt) F d(j) d[[3 — tfl'pa]), 

(4.4) 

where the perturbative Euler potentials a,/3 are now in¬ 
variant under a time-redefinition t ^ t + C. Under this 
choice the force-free equations result, 

da A d* Fb + dtp Ad* 6F = 0 , ( 4 - 5 ) 

{d/3 — ail'pdt) Ad* Fb 

+ {dip2 — ^pdt + d(/>) Ad* SF = 0 . ( 4 - 6 ) 

which contain no explicit t dependence. 

Note that these equations are coupled and non- 
separable. To proceed our analysis, we next adopt some 
simplifying approximations which enable us to handle 
them. These approximations will allow us to study their 
modal structure, under different limits, in the next sec¬ 
tions. 


B. Eikonal limit perturbations in flat spacetime 


We begin with the simple case of a flat background 
spacetime. Despite its simplicity, it is already a good ap¬ 
proximation when we examine that part of the jet that 
lies far away from its host black hole. In addition, we 
shall leave the functions "^2 and idp generic, without 
specifying their detailed forms, but we do assume that 
the flux function ip is just /2, which is consistent with 
the class of solutions described in Sec. El If one wishes 
to consider the jet of, say, a rapidly-spinning black hole 
in the far-zone in the future, it will be straightforward to 
insert the corresponding ip{p) and repeat the analysis in 
this section. 

Note that as the perturbations are in general not sta¬ 
tionary and/or axi-symmetric, we do not have simplifica¬ 
tions to the perturbing Euler potentials like those leading 
to Eq. (2.41. Nevertheless, as the background jet solu¬ 


tion is axis-symmetric, we can Fourier expand a and p3 
in Eq. (|4.4|) in cp\ 


= E- 


^am{t,r,e), I3 = Y,F'^^ 


i{t,r, 0 ), 


(4.7) 

and solve for a^n, Pm for individual m’s separately with¬ 
out worrying about the background introducing coupling 
between different m’s. The resulting wave equations (c.f. 
Eqs. (4.51 and (4.6)) are still coupled (between a and /3) 
and non-separable, with the following form 


The quantities Hi 2 and Vi _2 are lengthy differential op¬ 
erators which we do not present until further restrictions 
are introduced. In order to make further progress and ob¬ 
tain physical solutions, we focus on the short-wavelength 
(compared to the size of the black hole) perturbations, 
which are relevant for most astrophysical scenarios. Un¬ 
der this eikonal limit, one writes am. Pm as 


iS/e 


Br, 


o^S/e 


(4.9) 


Here Am and Bm are the amplitude functions which are 
assumed to vary slowly with location (on the length scale 
of the black hole size), while S is the common phase 
factor shared by both wave components, and it changes 
much faster than Am and Bm under the eikonal approx¬ 
imation. To keep track of this separation of scales, we 
introduced the book-keeping symbol e which in essence 
labels the WKB orders. Similar WKB treatments for 
uncoupled waves can be found in Refs. [43l [66U68] . and 
readers interested in further details can consult these ref¬ 
erences. Note that the phase matching between am and 
Pm follows from the assumption that we are solving for a 
single eigen-mode, in which case phase coherence should 
be preserved during propagation. In general though, the 
field perturbations could comprise a linear combination 


of a multitude of different modes, where Eq. ((4.9)) is of 
course no longer directly applicable. 


With Eq. (4.9) plugged into Eq. (4.8) and only keeping 


the leading order terms in e, the coupled wave equations 
become a matrix equation of the form 


'Hi vr 


Am 

_V2 H2_ 


Bm 


= 0 , 


(4.10) 


where the detailed expressions for "Hi, 'H 2 , Vi and V 2 
are given in Appendix]^ The determinant of the above 
matrix has to be zero for there to be a solution for Am 
and Bm, which requires that 'H 1 H 2 = V 1 V 2 (henceforth 
refereed to as the “determinant equation”). After some 
lengthy but nevertheless straightforward manipulations, 
we can factorize the determinant equation into the prod¬ 
uct of the following Hamilton-Jacobi equations 


^±Ha,b = Q, 
at 


(4.11) 


with 




Hu = 


47r2p2 J 


-np/ip) + m]± 




where the “momenta” are defined by 

_ dS 

PC = 


, (4.12) 


Hiam + ViPm — 0 , 
H2Pm + = 0 . 


(4.8) 


(4.13) 





















with C being any coordinate. We also note that = m 
and the terms inside the square bracket in Hb is pro¬ 
portional to — E^, which is a non-negative geometric 
scalar for the force-free plasma. These two Hamilton- 
Jacobi equations describe two different sets of modes. 
The first one is equivalent to 


= v^'^d^Sd^S = 0 , 


(4.14) 


which is the eikonal limit of the wave equation for vac¬ 
uum electromagnetic quasinormal modes. As discussed 
in Ref. |32j by writing down the local dispersion re¬ 
lation, this family of modes corresponds to the fast- 
magnetosonic waves in the force-free plasma; it is also 
consistent with the (eikonlal limit of the) “trapped mode” 
classification we found in analyzing the modal structure 
of the Blandford-Znajek solution m- We emphasize that 
this agreement is only in the eikonal limit, as the trapped 
modes differ from the vacuum electromagnetic quasinor¬ 
mal modes for longer wavelengths. 

The propagation of the second family of modes can be 
determined through the Hamiltonian equations of mo¬ 
tion. In fact, based on the form of Hb in Eq. (4.12), we 
can see that 

dz OHb 

dt dpz 


1 


d(j) OHb 
dt dp^ 

1 

(l + 4^) 




-^F T 7 ^\I^ + I > 


27rp^ y 


(4.15) 


where the ± sign labels the propagation direction and 
the equations above suggest that this wave has a null 
group velocity (as opposed to phase velocity) and freely 
propagates along the magnetic field lines in the frame of 
vanishing electric field (i.e. the “rest-frame” of the mag¬ 
netic fields). These properties are consistent with those 
of the Alfven waves in the force-free plasma [32] . On the 
other hand, this family of waves can obviously be seen 
as the eikonal limit of the “traveling waves” identihed in 
Ref. [37] (in that paper the travelling waves follow the 
mostly radial monopole-like background magnetic field 
lines). 

Now that we have introduced the necessary WKB tech¬ 
niques for handling coupled wave equations and learned 
about the mode structure of jets in the flat spacetime, we 


are now ready to study the wave propagation properties 
in a slowly-spinning Kerr background, where relativistic 
effects are non-negligible. 

C. Eikonal waves in a slowly-spinning Kerr 
spacetime 


With a rotating black hole sourcing the jet, we can 
still write down the wave equations of the jet perturba¬ 
tions based on Eqs. (4.5) and (4.6), subject to the Kerr 
background; 


ds^ = - 1- 


2Mr 


dH- 


AaMr sin^ 0 




dtdcj) 


+ ^dr"^ + p^dO'^ 
A 


-1- sin^ 9 I 


2Ma^r sin^ 9 


dcj )^. (4.16) 


Here A = — 2Mr -|- a? and p^ = cos^ 9. As 

discussed in Sec. [TT| the family of jet solutions is only 
known for a <C 1, in which case '0,12^;’,/ are expanded 


with respect to a and solved order by order (see Eq. 2.6). 


Under this slow rotation limit, we can still Fourier de¬ 


compose a and /3 as in Eq. (4.7), and the wave equations 
have the same form as Eq. (4.8), but with the operators 
now given by 


Hi = h: 


(0) 

i 

r(0) 


aH 


( 1 ) 




Oia^), 


(4.17) 


where 


= - 


dH 


f sin^ 9 


dcos9‘^ y dr 


= - 


r — M + M cos 29 + f ( ^ ^ 


rf 


df^ 


d9 


— cos9 — 
dr 


[(2/ — 1) cos 20 — 1] ^ -I- ^ sin20^ 


2r 


2r2 


50 


-^(cot2 0 + /) , 



2 f d 

= if sin 9m cot0— 
V 50 


im 

f n d 

cot 0——1 


sin^ 0 

V 50 

(notice that Hf''^ 

and 14^°^ 


+ fr 


dr 


fr 


dr 


(4.18) 


give the wave equations of 
the Wald background configuration, with a Schwarzschild 
black hole and surrounding force-free plasma), 
and 
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ijW = 






+ 2r2f}«-+2f^- 

dr dr ^ at de do 


+ I fr^ 




dj-2 


2/M- 


i9r 


+ / 


d'^'ip^^ 

ae^ 


- fcot{9) 


a^ 

90 


= 


^sin(0)|+2r2sin2(0) (frsm{9)^^+cos{e)^^\ ^ - 2/V sin2(0) cos(0)^^ — 


,a4^^ a2 


-2/Vsin3(0) 


a^ 

ar ^ 


+ 2/^rsin^(0) sin(0) 


ae 


- rcos(0) 


av’^^^ \ a^ 


-2/»„=(«) (2Af - r^n") ^ (,3s,m _ 


dr I drd9 
d^ 


2 sin(0) 


■ sin(0) 


a 

J at 


(4.19) 


^ / atar ' r v "“'v ata9 
-2M{l + 3cos{9))+4r^{r-3Msm‘^{9))n^p'' +r^{r-2Msi-a^{9))sm‘^{9)n^p^ 
sin(20)/(8M - 5r)^ + 2(3cos(20) + l)fr^^ + fr^smi29)^^^ - 2 fr sin^ {9)^"'^^ 


( 1 ) 


d9 


dr 


dra9 


-2sm^{9) f{5r-2M)cos{9)'^ -2/rsin(6»)^^^|-+/^r (rcos{9)^^-^ -sin(6l) 


( 1 )' 


+2/m^ f/rsin(6»)^^|^ 


+ cos(9) 


dip^'' 

a9 


a9‘^ 


d_ 

d9 


d^ 

dr 


(4.20) 


= - 


AimM{M cos(20) — M + r) a 


/^4 


at ’ 


= 


1 j 2r4 
2r6 1 / 






,3 ( aF2^’ 


/ 

Ar* 

+ — rM(cot^(0) — 1) + r^(3M — r csc^(0))l — 

/ at 


+r^ 


+2r^ 


4/r 


a2?)>F^ 


a2 


(4.21) 


, a# ' _ 

dr I drd9 


d9^ 


^ - 4cot(9)(5M - 2r)^^ - fr^ csc^(9) f 2 sm(29) ^ + (3cos(26») + 1 )^^^ 


2/r^ cot(0) 


a^V'F) 

dr"^ 


a9 


drd9 


dr 


dra9 


+2m^r^ csc^(0) /r 


dr 


cot(0) 


aV' 


( 1 )' 


ae 


ar 


a9 


a 


_a 

ar 


(4.22) 


While the equations are generally coupled in their a 
and a variables, and their r and 9 dependence are not sep¬ 
arable, it remains instructive to apply the eikonal approx¬ 
imation, and study the propagation of short-wavelength 
wave packets. 


Using the WKB techniqu e illu strated in Sec. |IVB[ We 
first make the ansatz of Eq. (4.9), obtain the matrix equa¬ 


tion as in Eq. ( |4.10[ ) but with different matrix compo¬ 
nents (see Appendix [B), and then solve the determinant 
equation. After some lengthy calculations we can show 
that two factor equations for dtS follow naturally from 
the determinant equation (up to 0{a)). The first one 
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FIG. 5: A demonstrative example showing the trajectories 
(orange curves) of the wave packets inside the jet region 
(bounded by the transparent tube), following the equation 
of motion as determined by Eq. | |4.24 L For this example, the 
mass of the black hole is chosen to be 1 and the spin 0.1. The 
value of flp is fixed at QhI^ and I is subsequently given by 
Eq. ( [2T0| ). 


reads 

/ dS\^ AamM dS 2 2 ^sc^ 0 fpg 

\ dt ’ 

(4.23) 

which is just the vacuum Teukolsky equation (before sep¬ 
arating out the 6 dependence) in the eikonal limit 
Therefore, clearly this family of modes are the “trapped” 
modes in the eikonal limit, as we have seen in the flat 
spacetime case. 

The second equation contains / and Vlp, so that this 
family of modes is affected by the background electro¬ 
magnetic field, in contrast to the trapped modes. The 
eikonal wave equation in this case reads 


Impz — r) — f^Mr'^ cos(20)] 

- sin^(6») ^2mflF + 

-b irfFr^ sin^ Opl = 0 , (4.24) 


with 


F = r — 2M sin^ 9, Pz = r cos 9pr — sin 6pg . (4.25) 


We note that Eq. (4.241 suggests that the wave prop¬ 
agates along z and (p directions. It is thus apparent 
that the second family of modes are traveling waves, as 


^ See equation 4.7 in m and keeping only the 0(a) terms 


expected. We present some example trajectories of wave 
packets following the equation of motion determined by 
Eq. ( frlit in Fig. 


D. Stability of the jet solutions 


As shown in Sec |III[ given an external magnetic field 
and a rotating black hole, the relaxed state of the mag¬ 
netosphere is independent of the initial configuration 
we chose. On the other hand, we have also shown in 
Sec|ll]that there is a whole family of stationary and axis- 
symmetric solutions which satisfy the force-free condi¬ 
tions. To reconcile the uniqueness of the steady state 
with the large available pool of candidates, it is natural 
to suppose that of all the members of the jet-solution 
family, only the steady state is in fact stable. 

To ascertain if this were the case, we need to exam¬ 
ine the stability of the family of solutions, with modal 
stability being the simplest approach. The analysis in 
Sec jlV C| and Sec IV D focuses on the wave equations in 
the eikonal limit, where wave damping or amplification 
may be ignored as they are slow when compared with 
wave oscillations. In order to study their stability, we 
need to go beyond the leading order WKB analysis em¬ 
ployed thus far. This is a technically arduous task to 
carry out, and beyond the scope of this work, but we can 
offer some physical arguments to this end. First, since we 
only impose axis-symmetric initial data in Sec |IIH and 
such symmetry is well preserved dynamically, it is suf¬ 
ficient to concentrate on the modes with m = 0. For 


these modes, Eqs. (4.17), (4.18) and (4.19), suggest that 


a “almost” describes the propagation of trapped modes 
while P “almost” describes the propagation of traveling 
waves. This observation follows from the fact that for 
the the Schwarzschild background (a = 0), a and /3 ex¬ 
actly decribe trapped and propagating modes. For the 
slowly-spinning-Kerr background, a and /3 fields are are 
only midly coupled with each other since the off-diagonal 
terms (given by Vi) are of order 0{a), while they prin¬ 
cipal operators Fli change only at order 0{a'^). Conse¬ 
quently the resulting eigen-frequencies can be at most 
order O(a^) away from their Schwarzschild limits |69) . 
Now, recall that trapped modes decay in Schwarzschild 
therefore trapped modes for small values of a can not sud- 
dently become unstable. On the other hand, traveling 
waves in the Schwarzschild spacetime should have zero 
damping as their scattering potential is infinitely shallow. 
Therefore, their counterparts in a slowly-spinning-Kerr 
spacetime can be unstable with a growth rate ^ 0{a'^) at 
most. However, this timescale is apparently much longer 
than the transient period we observe in Sec |III[ The argu¬ 
ments above thus point to modes, if unstable, only being 
so in much longer timescales than those observed in the 
dynamical studies. 

Of course, we emphasize that the above arguments 
are based on modal stability, which is different from lin- 
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ear stability. While our numerical investigation disfavors 
modal instability as the feature that renders alternative 
steady state candidates unfavorable, it does not also rule 
out non-trivial linear instabilities. However, given the 
broadness of the range of available jet-solutions, our ex¬ 
pectation is that this “selection rule” comes instead from 
more general mechanisms, such as the physics of current 
sheets -places where the force-free equations break down 
as they signal electrically dominated regions. 


V. DISCUSSION 

We have presented jet-like solutions of the force-free 
equations for spacetimes containing a slowly-spinning 
black hole. Our discussion illustrates that a family of 
such solutions can be found. Interestingly though, our 
numerical simulations clearly indicate that the configu¬ 
rations tend to change in time and approach a unique fi¬ 
nal steady solution, regardless of the initial configuration 
chosen. We have also carried out further analysis that 
disfavours modal instability as the triggering mechanism 
for this dynamical development. Instead, such behavior 
appears connected to the development of an equatorial 
current sheet which triggers Alfven waves propagating 
through the tube and re-arranging the solution towards 
the final one (similar observations have been made in the 
context of neutron stars [171I7SC5]). 

This would imply that current sheets play a vital role 
in the dynamical evolution of force-free magnetospheres. 
It dissipates a portion of the field energy, dynamically 
feeds back into the field evolution and helps to stabi¬ 
lize the magnetosphere. Therefore a complete study of 
the stability of different magnetospheric configurations 
should include current sheets as dynamical entities, in 
addition to the force-free plasma. An analytical treat¬ 
ment of this kind is still lacking, which will be the subject 
of our future studies. 

We note that one striking consequence of this 
numerically-observed dynamical behavior is that the 


force-free solution in black hole spacetimes might also 
obey a sort of “no-hair theorem” whereby the final sta¬ 
ble (stationary) solution is uniquely determined by the 
external field, the black hole mass and its spin. This is 
natural to be expected for Einstein-|-Maxwell black holes, 
but rather nontrivial with the presence of plasma and 
current sheets as we discussed. 

As a final comment We note that quite differently 
from tracing the dynamics of individual plasma particles, 
which is usually the topic of study for magnetospheric 
radiations, the “modes” investigated in this work repre¬ 
sents phonon-like collective motions of the entire mag¬ 
netospheric plasma. Such collective modes encode the 
fundamental parameters of the black hole and the large 
scale structures within the magnetosphere and may be 
tied to potentially observable phenomena. 
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Appendix A: The flat spacetime force-free equations 
in the eikonal limit 


The operators appearing in Eq. (4.10) are 
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2 r^Pr^^ ( + iticsc^{ 9 ) ) + ‘ 2 pe-^ ( + mcsc^{ 6 ) ) + 2 mr^VtpCSC^{ 9 ) 


dS d4’2 


dt 


89 


d'ip2 


89 


8S 


8t 


+rVr (^-(^^y+rV„-csc^m'j -pI ('r* - r^nj. + csc^(e)' 

(^7) +(w) -'■’‘*1)}' 


•Ha =-r2 -r2sin2(6») 


95 ; 

lit 


+ ml + sin^(0)pg + sin^(0) cos^{ 9 )p^ — 2r sin^(0) cos{ 9 )pgp. 


). 




Vi =rsiii(0)pr ^r^0i;’Sin^(6») + sin(0)pe ^sin(0)^^ +rcos(0)^^^ + m 


8S\ 


+ cos{ 9 )pg ( r^flp sin { 9 ) ( -^ ) + w ) + ( rsin( 0 )-^ + cos{ 9 ) 


8t J 


,^^rps(9i^'\ (r^sin^(e)f^y- 


89 J 


8t 


- sin'^{ 9 ) cos{ 9 )pl^ - rsin^{ 9 )pl^ , 

Va =-rpr (^r‘^npsm^{9) + sin(6<)pe ^sin(6i)^ +rcos(6»)^^ + m 


-Pe \r 


I 2 • /n\ /a\ 2^V'2 . 2/n\ 2^1‘2 

+ r sin(6>) cos{ 9 )p^ + r sin ( 9 )pg . 


55' 

lit 


(Al) 

(A2) 


(A3) 


sm(6i) cos(6») +mcot{ 9 )^ + + cot(6')^^ sin^(6') ^ 


(A4) 


Note that we have defined the momenta to be 


PC = 


dC’ 


(A5) 


for any coordinate C. We also note that we have the 
relationship 


dil’2 


8lp2 


I = — 2 TTrsin‘^{ 9 ) ( sin(0)-^^-rcos(0)^— ) (A6) 


8 r 


from Eq. (2.7), which is utilized when deriving Eq. (4.121. 


Appendix B: The Kerr spacetime force-free 
equations in the eikonal limit 

The operators appearing in the Kerr spacetime version 
of Eq. (|4.10|) are 
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Hi = 


csc(0) 

/^2 


2 2 ^V’2 ^ 


—2afmrpr 


dr 


f) f ^'^2 ^ I r, rl(l)^ I f2 2 2 I 1- 2'\ 

2afmpe—^ - r + ^amViy \ + f r p^ + fpg\ , 




(r — 2M sixr{9)) \—AamM s\rr{9)—— r'^ sur{9) { —] +fm^r\ -\-f^r'^ siic{9)pg 


+/^r^ sin^(20)p^ — 8/^r^sin^(0) cos{9)pgpr | 
Vi =m {frsm{9)pr + cos{9)pg) 


a I dip. 


( 1 ) 


2fr d9 


2cos(0) fm^ — r'^ siii"^(0) 


as; 

1h 


+ f‘^r^ sin(20) sin(0)p^ — 2f^rsiir{9)pgpr 


!sin(6*) ^2M — (fr sin{9)pr + cos{9)pg) 


'm 


- sin(0) 

^ cos(20)) ( ^ ) + /r^ sin(20)pgpr — 2frs\T?{9)pg + Arn^M — 2m‘^r 


( 1 ) 


dr 


dS 


dt 


dS 


dt 


V 2 = - m (/r sin(6>)pr + cos{9)pg) 


a I dip. 


( 1 ) 


2fr d9 


dS 


dt 


2cos(6>) fm^ — r^ sin^(0) — ) + f^r^ sin(20) sin{9)p^ — 2f^r sm'^{9)pepr 


+ sin(0) 
-fr 


2sin(6i) { 2 M - {frsin{9)pr + cos{9)pg) ^ 


dip2'’ / 3 


d \ I f ^ ) + fr'^ sm{29)pgpr — 2frsm'^{9)pg + Am'^M — 2Tn^r 


(Bl) 


(B2) 


(B3) 


(B4) 


where we have expanded the expressions up to 0{a), their terms cancel out in the determinant equation, which 
adopting expanded expressions such as 17^^ for clarity. written as 

Despite the unwieldy appearance of Hi and V^, many of 


dSV 


^7rr sin2(6») (j^pl - r^ (2amn^P + ~ csc’^i^) ) - fr {fpl + q^) 


—AamM csc^{9) 
,dS 


- r)p, ^sin2(6l) ^/r (fpl + qfj - + /m^F^ 


+7r/V2 sm\9)plF - — ) sin4(0)F2 


dS 


(B5) 


r 


It is now straight-forward to verify that the product of equation above at 0{a). 
Eqs. (4.23), (4.24) and agrees with the determinant 
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